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Abstract 

We show that oscillations are excited in a complex system under the influence of the external force, if the 
parameters of the system experience rapid change due to the changes in its internal structure. This excitation 
is collision-like and does not require any phase coherence or periodicity. The change of the internal structure 
may be achieved by other means which may require much lower energy expenses. The mechanism suggests 
control over switching oscillations on and off and may be of practical use. 

PACS numbers: 45.05,+x, 45.30,+s, 89.90.+n 
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A complex system subjected to an external load is a quite common phenomenon in nature, 
laboratory, and industry. Almost all systems experiencing the influence of external force possess 
internal structure which is not affected by this force directly but may change because of intrinsic 
dynamics of other external causes. The examples of such systems span all possible scales. A 
planetary magnetosphere under the influence of the solar wind is in the equilibrium (defined by 
the stationarity of the magnetopause position) if the time-independent incident plasma pressure 
is balanced by the magnetic pressure [|l]j. Yet the "stiffness" of the compressed magnetosphere 
is determined, among others,by the magnetosphere-ionosphere-atmosphere-solid earth interaction 
[@] . A muscle can be heuristically represented as a set of springs connected in parallel and in series, 
the elasticity of these "springs" being dependent on chemical processes and electrical signals in 
muscular fibers [TJ]. 

These systems may oscillate globally (the position of the magnetopause, the length of the mus- 
cle) near the equilibrium position (if not overdamped). There is a rapid growth of interest in 
possible effect of the internal variation on the excitation of global modes. Tanimoto and Um 
|Q] suggested that interaction with dynamical atmosphere and releases of energy due to weak lo- 
cal earthquakes may cause persisting global earth oscillations. Kepko and Kivelson [|5|] find a 
relation between magneto spheric bursty bulk flows and low frequency Earth magnetic field oscil- 
lations. Recently, Chagelishvili et al. ^ proposed a mechanism of oscillation excitation in a one- 
dimensional oscillatory system, based on the rapid change of the eigenfrequency. Unfortunately, 
their analytical consideration was inappropriate which obscured the physics of the phenomenon 
and lead them to conclusions which are too restricted. 

In the present paper we investigate the excitation of oscillations in complex systems due to 
rapid changes in internal parameters under quite general conditions. We show that such excitation 
(and oscillation amplification) is a quite common phenomenon and calculate its efficiency. This 
question may be of significant importance in applications. One can easily imagine systems where 
such excitation is unwanted (like building constructions). On the other hand, the possibility of 
generating oscillations by not changing external force but applying only weak forces (and energy) 
to cause local rapid changes of the system parameters is rather attractive. A simple system where 
both situations are possible is an LC circuit with distributed capacity and connected to a constant 
emf. Excitation of oscillations at the circuit frequency is unwanted when it is used in a device 
designed to measure electromagnetic spectra. On the other hand, such excitation would be useful 
in sustaining certain level of the current when weak damping is present and it is difficult to apply 
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periodic forcing. 

In what follows we consider the systems described by the following (vector) equation of motion 

dU 

MX = -— + F, (1) 

where X is the vector of external parameters of the system (position) which are subject to the 
external force (position of the magnetopause under the solar wind pressure, length of the loaded 
muscle, charge on the effective capacitor when the LC circuit is connected to a constant emf). 
Let also x be the vector of internal parameters which are affected in a different way or vary at 
a much small time scale for some reasons. The generalized mass matrix M and the potential U 
depend both on the external and X and internal x coordinates, and the external force F depends on 
time. When the equations governing the behavior of the internal coordinates x are not known or 
are too complicated we may phenomenologically describe their influence as temporal dependence 
of M and U. If the system parameters and external force do not vary with time, the system is 
in the equilibrium position X e(? which is determined by the condition dU / <9X eg = F. Near the 
equilibrium the potential can be written as 

C/=iK-(X-X m )-(X-X m ), (2) 

where K is the stiffness matrix, so that one has 

MX = -K(X - X m ) + F, (3) 

which looks exactly as a usual oscillator equation except that now mass and stiffness are matrices. 
Recently it was shown that the process of vortices generation in magnetohydrodynamic and shear 
flows [0] can be described by a special one-dimensional case of Eq. (Q). 

The equilibrium position is X eg = X m + K _1 F, and the general solution of Eq. (Q) near the 
equilibrium point is 

X = X eq + 2j a%h sm(uit + (pi), (4) 

i 

where the frequencies tOi and unity vectors are the eigenvalues and eigenvectors of the matrix 
W = M~ l K, respectively. Since the equilibrium is assumed to be stable, the matrix W is posi- 
tively determined and all frequencies are real. a 4 and 0, are the amplitudes and initial phases of 
the normal modes. 
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Mostly known channels of the energy input to the oscillating system include adiabatic change 
of the natural frequency, resonance with the periodically changing external force, or parametric 
resonance (see, for example, Ref. [|8p). The last two imply (quasi)periodic behavior of the nat- 
ural frequency or external force and are not within the scope of the present paper where only 
nonperiodic changes are considered. 

The motion of the system is now fully determined by M, K, X e(? , and F. If these parameters 
change slowly at the typical time scale of oscillations, that is, the typical time of variation is much 
larger than all 1/ui, the equilibrium adiabatically shifts its position, while the amplitudes follow 
the well-known adiabatic law ujiaf ~ const. In particular, the system which was in the equilibrium 
in the beginning will remain in the equilibrium in the course of the parameter variation, so that no 
oscillations are excited. 

In the present paper we consider the case of rapid parameter changes, where the typical time 
of variation T <C l/cu,, for all i. We show that there this, in general, results in the excita- 
tion/amplification of oscillations in the system. A numerical example of a special type of exci- 
tation in the simplest one-dimensional system was recently considered by Chagelishvili et al. [§]. 
Here we consider a most general case of nonadiabatic excitation of oscillations in loaded systems 
near the equilibrium point. 

To study quantitatively the effect let us assume that all parameters vary only within the time 
interval < t < T, and lu-T < 1. It is convenient to define W = M^K and f = M _1 F. From 
Eq. (Q) one immediately has 

X(T) = [ T [f (t) - W(t) (X(t) - X m (t))] , (5) 
Jo 

where the time dependence of the system parameters is explicitly shown. In what follows we 
denote the system parameters and variables at t = and at t = T with subscripts 1 and 2, 
respectively. The solutions at t < and t >T will take the form: 

Xfc(t) = X fce g + ^ a ki&ki Sm(u ki t + <f>ki), (6) 
i 

where k = 1 for t < and k = 2 for t > T. Since uj{T <^ 1 one can neglect the difference 
between X^T) and X 2 (T), so that one has 

i i 
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The second matching condition is obtained from Eq. (Q) in the following form: 
X 2 (T)-X 1 (T)= / (f-fi)dt 



^0 



(iy - w 1 )dt\x 1 (T) - [ (wx m - w 1 x lm )dt. 



(8) 



It is convenient to define the instantaneous equilibrium position as X eq = X m + W x f . Taking 
into account Eq. (Eh one has 



22 UJ 2ia 2 ie 2 i cos(uj 2 iT + </> 2i ) = 2j ^li^eii cos(co> H T + H ) 

+ [/ (W^-^OtftlVaiiew^wiiT + ^ + l/ W(X leq -X eq )dt}. 
Jo „• ./o 



(9) 



Eqs. ([7]) and (Q) allow one to find a 2i and 2?; knowing the initial state. 

It is easily seen that even when an = for all i, the oscillation amplitude in the final state is, 
in general, nonzero: 

4 = [^(X 2e(? - X leq )} 2 + -![ f e* 2l W{X - X leq )}\ (10) 

Thus, the oscillation is excited due to the (a) irreversible shift of the equilibrium position (first 
term in Eq. (|T0|)) and (b) reversible temporal shift of the equilibrium position with subsequent 
return to the same equilibrium (the second term). The nature of the energy input is different for 
the two mechanisms (instantaneous change of the potential energy in the first case, and work 
done by the external force because of the excursion of the system from the equilibrium in the 
second case) but in both cases the interaction is collision-like: certain amount of momentum and 
energy is transferred to the system in a short time interval. The physical nature of the excitation is 
quite different from the proposed earlier [Q] quasi-parametric resonant interaction with a Fourier- 
component of the changing eigenfrequency. It is worth noting that although ( |TO| ) depends on the 
change of all parameters, M, K, and X m , the change of the mass M along does not affect the 
instantaneous equilibrium position X eq and therefore does not result in the oscillation excitation, 
as could be expected. 

Within the chosen approximation Eq. dlOt ) includes all possible internal perturbations of M, K, 
and X m , and external perturbations of F, thus giving the most general description for nonadiabatic 
excitation of complex systems. It includes the system considered in Ref. |] as a special one- 
dimensional case where only K is changed. It is instructive to rewrite Eq. ( JT0| ) for this case, where 
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X m = 0, and fi = f2,W = u) 2 , and u\ = u>2- Then the excited amplitude takes the following 
simple form: 



a 2 



T cu 2 



fi f 



<jj- 



dt 



(ID 



l UJl \UJ± 

from which one can easily see that the excitation occurs always and not only when the eigenfre- 
quency decreases in the perturbation (cf. Ref. ||). 

In what follows for simplicity of presentation we restrict ourselves with the one-dimensional 
case. Multidimensional generalization is straightforward. In the one-dimensional case Eqs. (j7|) 
and (SI) immediately give 



a\ = [{X leq - X 2eq ) + ai sin(cJiT + 0i)] 2 
+ 0J2 2 [aiUJi cos(dJ 1 T + (j)i) + ( / (oj 2 — uj\)dt)ai sin^T + 



(12) 



+ / co 2 (X leq -X eq )dt} 2 . 
Jo 



In most cases the phase <f>x of the interaction beginning is unknown (unless the perturbation is 
carefully prepared with some definite purpose in mind). In this case one can consider the phase </>i 
as random and average over random distribution to obtain eventually: 



2 U 1 



1 uj 2 qV-^K) 



(13) 



leg 



Xop.a) H ^ 



^ 2 (Xl eq — X eq )dt 



Eq. ( JT3| ) describes the amplification (reduction) of oscillations (first term) and excitation of os- 
cillations due to the shift of the of the equilibrium position and to the work of the external force 
during temporal shift of the equilibrium. 

In the particular case of perturbation where after time interval T the system parameters return 
to their initial values, one arrives again at 



(Xc\ CL-\ 



T u 2 



J U)\ \u 



-f4--4u 



(14) 



where we neglected the term containing [J q T (uj 2 — uj\)dt} 2 /uo 2 <^ 1. Roughly speaking, during 
the excursion of the internal parameters from their equilibrium values the energy of oscillations 
increases by 5E = [J Q T uj 2 {X\ eq — X eq )dt} 2 on average. If the system experiences a series of iV 
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such rapid variations with randomly distributed phases, the total oscillation energy increase would 
be about N5E, without any necessity to arrange the phases or periodicity of the variations. 

Another efficient method of excitation is the rapid shift from the equilibrium position with 
subsequent return to this position after time ^> 1/uj, In this case, neglecting the last term in 
Eq. (JBj) and after some algebra one obtains 



which shows that this scenario always results in the oscillation amplification. 

To conclude, we have shown that nonadiabatic changes in the system parameters and/or exter- 
nal force are efficient in excitation or amplification of oscillations in driven oscillatory systems 
under external load. Randomly distributed in time, short nonadiabatic pulses result in efficient 
transfer of energy to the system. The energy transfer manifests itself in continuously increasing 
oscillation amplitude. This amplitude increase is not restricted to the periodically repeated (co- 
herent) perturbations, as was suggested in [§J, but occurs in quite general conditions. The energy 
input effect is essentially nonresonant and more collision-like where additional momentum/energy 
are transferred to the system at the time scale much smaller than the typical timescale of varia- 
tions in the system. The effect may be important for the systems whose natural frequencies may 
vary quickly due to the variable internal coupling. It should be emphasized that the importance 
of the above analysis is well beyond the consideration of simple oscillatory systems, desribed by 
the simple oscillatory potential in the form (0) (chosen here for convenience and simplicity of 
presentation), but may be applied to quite general systems. The results (qualitatively) are valid for 
any system capable of (generally nonlinear) oscillations near the forced equilibrium position (al- 
though quantitative calculations would require knowledge of the structure of a particular system 
and ability to translate it into a low-dimensional description with small number of parameters). 
There is a wide spectrum of such systems, from large astrophysical scales (gravitationally bound 
systems, planetary magnetospheres under solar wind influence) to usual human scale (muscles, 
constructions) and down to small scales (electric circuits). Such generation of oscillations may be 
unwanted in some systems, like possible excitation of internal currents in spacecraft circuits by 
cosmic rays. On the other hand, the essentially nonresonant generation of oscillations might be 
useful in experimental determination of the natural frequencies of the systems where it is difficult 
to apply periodic external force but where the internal structure can be changed relatively easily 
and rapidly. Such methods can be also used to re-excite damped oscillations without changing of 




(15) 
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the main load. 

Finally, let us use a very simple model to see whether reconnection at the dayside magnetopause 
may be responsible for excitations of global magneto spheric oscillations. The position of the 
magnetopause is determined by the balance of the incident plasma pressure n u m p V 2 (where n u 
and V u are the solar wind density and velocity, respectively, and m p is the proton mass) and the 
magnetic pressure B 2 /8tc. Is the magnetopause is compressed by x, the magnetohydrodynamic 
conservation of magnetic flux [§] predicts that the magnetic field increases as L/(L — x) where L 
is the equilibrium standoff distance of the magnetopause. Thus, there appears the excess force of 
~ 2n u m p V 2 xA/ L, where A is the effective area. This force has to accelerate the mass of about 
n d m p AD, where n d » 5n u is the average plasma density, and D is the distance between the shock 
and magnetopause. Using the typical parameters V u ~ 400 km/s, L ~ IORe, and D ~ R E , 
where Re ~ 6, 000 km is the Earth radius [Q], one finds the typical oscillation periods of the order 
of ~ 10 min. On the other hand, the typical time of reconnection should be of the order d/v n , 



where d ~ 800 km is the magnetopause width [|10|], while the velocity v n may be as high as 50 
km/s [[TT|]. Thus, the typical time of reconnection ~ 10 sec and much smaller than the oscillation 
period. Reconnection results in the breakdown of magnetohydrodynamic s and therefore reduces 
the "stiffness" of the magnetic field, thus effectively temporarily reducing the global oscillation 
frequency. Hence, the conditions of Eq. (|T^) are satisfied and excitation is possible. Of course, 
quantitative calculations require that we are able to translate the reconnection process into the 
change of internal parameters, so that at this stage the proposed scenario should be considered as 
a speculative hypothesis. 
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